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Abstract 



In this paper we introduce a new lifetime distribution by compounding ex- 
ponential and Poisson-Lindley distributions, named exponential Poisson-Lindley 
distribution. Several properties are derived, such as density, failure rate, mean 
lifetime, moments, order statistics and Renyi entropy. Furthermore, estimation by 
maximum likelihood and inference for large sample are discussed. The paper is 
motivated by two applications to real data sets and we hope that this model be 
able to attract wider applicability in survival and reliability. 
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2 1 Introduction 



Adamidis and Loukas (1998) introduced a two-parameter lifetime distribution with de- 
creasing failure rate by compounding exponential and geometric distributions, which 
was named exponential geometric (EG) distribution. In the same way, Kus (2007) and 
Tahmasbi and Rezaei (2008) introduced the exponential Poisson (EP) and exponential 
logarithmic distributions, respectively. 

Recently, Chahkandi and Ganjali (2009) introduced a class of distributions, named 
exponential power series (EPS) distributions, by compounding exponential and power 
series distributions, where compounding procedure follows the same way that was pre- 
viously carried out by Adamidis and Loukas (1998); this class contains the distributions 
mentioned before. Extensions of the EG distribution was given by Adamidis et al. (2005), 
Silva et al. (2010) and Barreto-Souza et al. (2010), where the last was obtained by com- 
pounding Weibull and geometric distributions. A three-parameter extension of the EP 



distribution was obtained by Barreto-Souza and Cribari-Neto (2009). 

Sankaran (1970) introduced the discrete Poisson-Lindley distribution by mixing Pois- 
son and Lindley (see Lindley, 1958) distributions. Ghitany et al. (2008) studied the 
zero-truncated Poisson-Lindley distribution and its applications. A random variable N 
with this distribution has probability function given by 

where 9 > 0. In this paper, we introduce a new lifetime distribution by compounding 
exponential and Poisson-Lindley distributions as follows: Let Zi, Z 2 , . . . , Zj^ be a random 
sample following an exponential distribution with scale parameter (3 > and probability 
density function in the form g(x) = j3e~ l3x , for x > 0, where N is a discrete random 
variable following Poisson-Lindley distribution with probability function given by (flj. 
The marginal cumulative probability function of X = min{Z i }^ 1 reduces to 

9 2 e~ fSx 1 + 6 + (2 + 6) (1 + 9- e-P*) 

F(x) = l-- — — — i ^ = '- x>0, 2 

K J 1 + 39 + 9 2 {i + e- e -^f 

which defines the exponential Poisson-Lindley distribution. We denote a random variable 
X with cdf given by (J2J) as X ~ EPL(/3, 6). We notice that, as Poisson Lindley distribution 
is not contained in class power series of distributions, our distribution is not contained in 
the class EPS of distributions introduced by Chahkandi and Ganjali (2009). To generate 
a random variable X ~ EPL(/3, 6), we computate X = i^f/) -1 , with U ~ U(0, 1), where 



F(^-^i og r (1+ ^" 1)+(6+ 

' '-^ 2 + 6 + a(l-U) 



-1//3 



with A v = [b-2a(l + 6)(U-l)] 2 + Aa(l + 6) 2 (U-l)[2 + 6 + a(l-U)},a = 6- 2 (l + 36 + d 2 ) 
and b = 3 + A9 + 9 2 . A simple interpretation of the proposed model comes from a situ- 
ation where failure (of a device for example) occurs due to the presence of an unknown 
number, say N, of initial defects of the same kind. The random variables Z's represent 
their lifetimes and each defect can be detected only after causing failure, in which case it 
is repaired perfectly. Thus, the distributional assumptions given earlier lead to the EPL 
distribution for modeling the time of the first failure. 

The paper is organized as follows. In Section 2, we derive properties of the new lifetime 
distribution, such as density, failure rate and mean residual lifetime. Expressions for the 
moment generating function and moments of the EPL distribution are presented in the 
Section 3. Further, order statistics and its moments also are discussed in this Section. An 
expression for Renyi entropy is given in the Section 4. Estimation by maximum likelihood 
and inference for a random sample from EPL distribution are discussed in the Section 5, 
while two applications to the real data sets are presented in the Section 6. Finally, we 
conclude the paper in the Section 7. 



2 Density function and failure rate 

The probability density and survival functions associated to (pi) are given by 

_ f36\l + 6) 2 e^ x 3 + 6 - e-? x 
f{X) ~ 1 + 36 + 6* (i + e-e^) 3 {) 



and 



6 2 e~P x 1 + 6 + (2 + 6) (1 + 6- e~ Px ) 
1 + 36 + 6 2 (1 + 6- e~^) 



S(x) = v^TT^ 1 , „ M fl _, 2 L , (4) 



for x > 0, respectively. We observe that pdf ^ can be expressed as 

/(*) = [ °l% \ g 2 \ffi(x) + 3g 2 (x) + g 3 (x)], 

where gi(x) = {(1 + 6)e l3x — 1}~\ for x > and i = 1,2,3. Of this way, density /(•) 
is a linear combination of monotone decreasing functions with positive coefficients and, 
therefore, it is monotone decreasing. Moreover, we have that linL r _> + f(x) = (3(8 3 + 
Ad 2 + 56 + 2)/(6 3 + 3d 2 + 6) and lim^^ f(x) = 0. We also have that f(x) -» /3e"^, 
when 6 — > oo. With this, the exponential distribution is obtained as limiting distribution. 
From (J3l) and Q, we obtain the failure rate of the EPL distribution: 

., . i3(i + e) 2 (3 + e-e-?')(i + e-e-'')-' 

h(x)= l + tf + (2 + tf) (!+«-«-») ' X> °- (5) 

The failure rate of the EPL distribution is decreasing according to the note of Proschan 
(1963) who proved that the decreasing failure rate property is inherent to mixtures of 
distributions with constant failure rate. Further, we have that lim a ._ > .o+ h(x) = (3(6 3 + 
Ad 2 + 56 + 2)/(6 3 + 36 2 + 6) and lim^^ h(x) = /3. Figure [I] shows plots of the pdf @ 
and failure rate (J5J) with (3 = 1 and some values of 6. 

Given survival to time xq, the residual life is the period from time xq until the time 
of failure. The mean residual lifetime of a continuous distribution with survival function 
S(x) is given by 



m(x ) = E(X — x \X > xq) = —. — r / S(t)dt. 

S{xo) J xo 

By taking S(x) given by Q, we obtain the mean residual lifetime of the EPL distri- 
bution: 

1 - (2 + 0)(1 + 6 - e~ px °) log[l - (1 - 6)- l e-P XQ ]eP Xo 

m ^> ~ f3[2 + 6 + {l + 6){l + d-e-^)- 1 } ' 

Based on the behavior of the failure function h(x) of the EPL distribution, which is 
decreasing, we conclude that m{x) is increasing (see Watson and Wells, 1961). 





Figure 1: Plots of the pdf and failure rate of the EPL distribution for some values of the 
parameters. 

3 Moments and order statistics 



As with any other distribution, many of the interesting characteristics and features of 
the EPL distribution can be studied through the moments. Let X be a random variable 
following EPL distribution with parameters /3 > and 9 > 0. We now calculate the 
moment generation function of X denoted by Mx(t) = E(e~ tx ). For \z\ < 1 and k > 0, 
we have that 



;i-*)-* = f; r !*+:?V 



3=0 



T(k)j\ 



(6) 



Using this expansion in the term {! — (! + 9) l e l3x } 3 and for t > —j3, it follows that 



M x (t) 



:i + 



\ — 1 /"OO 



i + 36 + e 2 



-{P+t)x_ 



3 + 9 - e~P x 



dx 



1 + 39 + 9 2 



{1 - (1 + 0)-l e -*} 



x 



and then 

M x (t) 



3=0 






£ 



(7) 



for £ > — f3. With formula (J7|), we can obtain an expression for the rth moment by 
calculating E(X r ) = (-l) r d r M x (t)/dt r \ t=0 . Hence, it follows that 



E{X r 



r\9 2 



2/3 r (l + 39 + 9 2 ' 



£(* + 



u+1) f (3 + 9)(j + 2) j + 1 



i=o 



(j + 1; 



(i + 2) r 



We can express the above formula in terms of the function L n (z) = Y^jLi zi lJ n i which is 
a generalization of Euler's dilogarithm function (see, Erdelyi et al., 1953, p. 31) and is 
also known as the polylogarithm function; the function is quickly evaluated and readily 
available in standard software such as MATHEMATICA and MAPLE. So, we obtain 

r\9 2 

E(X r ' 



v r _l 



1 



f3 r (l + 39 + 9 2 ) {"' x \9 
The mean and variance of X are given by 

9 



1 



(2 + 9)L r 



1 



+ 1 



E(X) 



/3(1 + 36» + 6» 2 ) 



l + 9(2 + 9)\og 1 + 







and 



Var(X) 



log (1 + |) 2 + 2(2 + 9)L 2 (^) 



1 + 39 + 9 2 



1 + 0(2 + 0) log(l + |) 
1 + 30 + 2 



1. 



respectively. Consequently, the coefficient of variation is 

Var(X) _ 2 (1 + 30 + 2 ) [log (l + |) + (2 + 9)L 2 (^)] 

EW 2 [l + 0(2 + 0)log(l + |)] 2 

Table [T] presents mean, variance and coefficient of variation for some values of the pa- 
rameters computed by using formula ((8|. We now discuss some properties of the order 



0| 


mean 


variance 


cv 


0.5 


0.4315028 


0.3923130 


2.107004 


1 


0.6158883 


0.5966291 


1.572898 


5 


0.9001530 


0.8989667 


1.109458 


10 


0.9494062 


0.9491146 


1.052966 



Table 1: Mean, variance and coefficient of variation for (3 = 1 and some values of 9. 



statistics from EPL distribution. Order statistics are among the most fundamental tools 
in non-parametric statistics and inference. They enter problems of estimation and hy- 
pothesis testing in a variety of ways. Let Xi, . . . ,X n be a random sample of the EPL 
distribution. The pdf of the ith order statistics X i:n is given by 



.//: 



X 



n\f3{l + 9) 2 9 2 ^ n 



-i+X) p -0x(n-i+l) 



(3 + 9 



7 -f3x\ 



\n 



i)\(i 



(1 + 9 



c\ln- 



1)!(1 + 39 + 9 2 ) n - i + 1 (l + 9- e -/te)2(n-0+3 
9 2 e -p x i + e + {2 + 0)(l + 9- e"^) 



[1 + 9 + (2 + 9) x 



1 + 39 + 9 2 



'1 + 



-0x\ 



i-\ 



for x > 0. We now derive an expression for the rth moment of the ith order statistic 
X i:n . For this, we use a result due to Barakat and Abdelkader (2004) and obtain 

k 



E{X\ 



E 



k=n—i+l 
oc 

J — 1 - 

x e 



i \k— n+i— 1 



9 2 (1 



\-i 



1 + 39 + 9 2 



-l e -/3a;|fc 



{1 



-Px\2k 



dx. 



Applying (pi) and binomial expansion to the terms {1 — (1 + 9)~ l e~^ x }~ and {3 + 9 — 
(2 + 9)(l + 9)~ 1 e~^ x } of the above equation, we obtain an expression for the rth moment 
of X.„: 



E(Xl 



T . 

w 

oo 

E 

3=0 



E E< 

k=n-i+l 1=0 



-1 



,fc— n+i+l— 1 



:i + 



n 



k\ 9 2k (3 + 6) k - l (2 + 6y 
l) (l + 9) l + k (l + 39 + 9 2 Y 



j + 2k-V 

2k -l )\k + i+ 3 y 



x 



(9) 



Table [2] lists some numerical values for the first four moments of order statistics X 1:2 o, 
Xlo:2o and -X" 2 o : 2o from equation ([9]), with the index j stopping at 100, and using numerical 
integration. The parameter values are taking as (3 = 0.1 and 9 = 0.5. The results in 
this table show good agreement between the two methods. We now are interested in the 



Xi-,20 4 r th moment—)- 



i = 1 Expression 
Numerica 



(9f 



i = 10 Expression 
Numerica 



(f 



% = 20 Expression 
Numerica 



(f 



0.1172126 0.03121316 0.01280193 0.007203924 
0.1260696 0.03281421 0.01323651 0.007361332 



2.279001 
2.056653 



4.940849 
4.899393 



13.43999 

13.42872 



42.12036 
42.11641 



22.10733 
22.10734 



621.2616 
621.2616 



21864.98 
21864.98 



944798.2 
944798.2 



Table 2: First four moments of some order statistics from ^ and via numerical integra- 
tion. 



asymptotic distributions of the extreme values X 1:n and X n:n . From equations ^ and 
Q, it can be seen using L'Hopital's rule that 

hm ' , = x and hm ^^^^ = e ' . 
t ^o+ F(t) «^oo S(t) 

From this and making use of Theorem 1.6.2 of Leadbetter et al. (1987), we obtain 



P (a n (X mn - b n ) <x)=e 



-/3x 



and P {cn{X Xm - d n ) < x) = 1 - e x , 



as n — > oo, where the norming constants a n ,b n ,c n > and d n can be determined by 
Corollary 1.6.3 in Leadbetter et al. (1987). For instance, we can select a n — 1 and 

b n = F(l — 1/n)^ 1 as n — > oo. 



4 Renyi entropy 

The entropy of a random variable X is a measure of uncertainty variation. The Renyi 
entropy is defined as Ir(o) = (1 — a)' 1 \og( J R f a (x)dx), a > and a^l. Let /(•) be 
the pdf of the EPL(/3, 9) distribution. Using expansion (pi) for (1 + 9 — e~ l3x ) 3a and after 



some algebra, we obtain that 

P a ~ 1 [9(l + 9)} 2a x 



F(x)dx- Ti3a)il + 3e + e2) 



E 

j=0 



T(3a + j)(3 



\2a+j 



X 



(1 + 6) 3a +ij\ 
1 — i?2+e (a, a + j — V 



3+0 



where B x (a,b) = f Q t a (l — t) b l dt is the incomplete beta function, for a, b > and 
x G (0, 1). Hence, an expression for Renyi entropy of the EPL distribution is given by 

I R {a) = -log/3 + (l-a)- 1 {2alog[fl(l + fl)]-logr(3a)-alog(l + 3# + # 2 )} + 



[l-a^log^ 

.3=0 



T(3a+j)(3 + 9) 
(1 + 9) 3a +ij\ 



2a+j 



1 — B2+1 (a, a + j — 1) 



3+S 



It is important to derive an expression for the Renyi entropy because it serves as a 
measure of the shape of a distribution and can be used to compare the tails and shapes 
of various frequently used densities; for more details, see Song (2001). 

5 Estimation and inference 

Let Xi, . . . ,X n be a random sample from EPL(/3, 9) distribution with observed values 
xi, . . . , x n . The log-likelihood function is given by 

£ n = £ n (f3,9) = n\ogj3 + 2nhge(l + e)-n\og(l + 3e + e 2 )-pJ2 Xi + 

i=\ 

•a n 

J2 log (3 + 6 - e-^) - 3 J^log (1 + 9 - e" fe ) . (10) 



»=i 



i=i 



The score function associated to log-likelihood (10) is U n = U n (j3, 9) = (d£ n /d/3, d£ n /d9) T , 
where 

1 3 



d/3~p 2^ x * + 2^ 



Xj,& 



i=l 



i=l 



3 + 9- e-fe 1 + 9 



-/3xi 



and 



d£ n 2n (1 + 29) n (3 + 29) 



09 



9(1 



1 + 39 + 9 2 



+ E 



i=l L 



3 + 9- e-fa 1 + 9 - e-to 



The maximum likelihood estimators of j3 and 9 are obtained by solving numerically 
the nonlinear system of equations U n = 0. It is usually more convenient, however, to use a 
nonlinear optimization algorithm (such as the quasi-Newton algorithm known as BFGS) 



to numerically maximize the log-likelihood function in (10). Moreover, under regularity 



conditions, the expected value of the score function vanishes. With this, we obtain 

1 3 \ 3 + 29 2 + A9 



E 



3 + 9- e~P x 1 + 9 



-fiX 



1 + 39 + 9 2 9 + 9 2 



and 



E 



Xe- 13 * 



3 + 6- e-? x 1 + 6 



-fix 



fl[l + 0(2 + 0) log (l + |)] 1 
f3(l + 36 + 6 2 ) /?' 



For inference, Fisher's information matrix is required. It is given by 

K n = K n (/3,6)=n 

where 

1 T X 2 e-? x 1 

3(1 + 6) E 



KfJ,/3 K fJ,e 



k^ = - + (3 + 6)E 



^e,e 



E 



(3 + 6-e-P x ) 2 
2{l + 26 + 26 2 ) 7 + Q6 + 26 2 
(1 + 30 + 2 ) 2 



X 2 e-? x 



6 2 (\ 
i 



+ E 



l + 6-e~P x ) 2 

2 



1 + 9- e-? x 



«e,/3 = E 



Xe-P x 



(3 + 6 -e 



3 + 6- e~P x 
-3E 



-fJX)2 



-3x 
Xe-ft* 



l + 6-e-P x ) 2 



The above expectations can be obtained numerically. Under conditions that are fulfilled 
for parameters in the interior of the parameter space but not on the boundary, we have 
that 



(5-(5 



An 



o ),*GM)~ ] ,- 



as n — > oo, where K((3,6) = lim. n _^ OQ n~ 1 K n (/3,6). The bivariate normal distribution 
with mean (0,0) T and covariance matrix i^ n (/3,0) _1 can be used to construct confidence 
intervals for the model parameters. 



6 Applications 

In this Section we fit exponential Poisson-Lindley distribution to two real data sets and 
compare with the exponential-geometric, exponential Poisson, exponential logarithmic, 
Weibull and gamma distributions, whose densities are given by 

fe g (x-(3 1 ,p 1 ) = (3,(1 - Pl ) e -^(l- Pl e-^)-\ ft > 0, Pl e (0, 1), 
f ep (x; fo, A) = _^_ e -A-/^+Acxp(-/^) ? £ 2> A > 0> 

1 6i(l - vo)e~^ 3X 

fei(x;P 3 ,p 2 ) = - P3{ P2) /3 3 >0, p 2 e(0,l), 

- logp 2 1 - (1 - p2)e-P- iX 

f w (x; /3 4 , a) = aP2x a ~ l e~ PAX i /3 4 , a > 0, 



fgfafai) 



r( 7 ) 



^ x i-i e -fcx^ /3 5;7 >o, 



for x > 0, respectively. The first data set is given by Linhart and Zucchini (1986), which 
represents the failure times of the air conditioning system of an airplane: 23, 261, 87, 7, 



120, 14, 62, 47, 225, 71, 246, 21, 42, 20, 5, 12, 120, 11, 3, 14, 71, 11, 14, 11, 16, 90, 1, 16, 
52, 95. 

In the second data set, we consider vinyl chloride data obtained from clean upgradient 
monitoring wells in mg/L; this data set was used for Bhaumik et al. (2009). The data 
are: 5.1, 1.2, 1.3, 0.6, 0.5, 2.4, 0.5, 1.1, 8.0, 0.8, 0.4, 0.6, 0.9, 0.4, 2.0, 0.5, 5.3, 3.2, 2.7, 
2.9, 2.5, 2.3, 1.0, 0.2, 0.1, 0.1, 1.8, 0.9, 2.0, 4.0, 6.8, 1.2, 0.4, 0.2. 

Tables |3] and [4] give us the fitted parameters, Kolmogorov-Smirnov (K-S) statistics 
and its respective p-values for the first and second data sets, respectively. We see that 
all the distributions in both tables show a good fit for the given data sets and, in each 
case, we do not reject the hypothesis that the data comes from distribution considered 
at any usual significance level. We also see that EPL distribution proves to be a good 
competitor for other distributions already known in literature. 



Distribution 


Estimates 


K-S statistic 


p- value 


EPL(M) 


(0.0101, 0.9193) 


0.1290 


0.6531 


EG(/3 1)Pl ) 


(0.0102, 0.6148) 


0.1262 


0.6793 


EP(/3 2 ,A) 


(0.0106, 1.7941) 


0.1472 


0.4890 


EL(/3 3 ,p 2 ) 


(0.0111, 0.1932) 


0.1288 


0.6555 


W(/9 4 ,a) 


(0.0183, 0.8533) 


0.1531 


0.4394 


G(/3 5 , 7 ) 


(0.0136, 0.8135) 


0.1694 


0.3187 



Table 3: Estimates of the parameters, Kolmogorov-Smirnov statistics and its respectively 
p-values for the first data set. 



Distribution 


Estimates 


K-S statistic 


p- value 


EPL(/3,0) 


(0.4796, 5.0811) 


0.0882 


0.9331 


EG(/3 ljPl ) 


(0.4818, 0.1771) 


0.0876 


0.9360 


EP(/3 2 ,A) 


(0.4767, 0.4276) 


0.0880 


0.9341 


EL(/3 3 ,p 2 ) 


(0.4867, 0.7022) 


0.0870 


0.9394 


W(/3 4 ,a) 


(0.5296, 1.0101) 


0.0918 


0.9116 


G(/3 5 ,7) 


(0.5654, 1.0626) 


0.0973 


0.8733 



Table 4: Estimates of the parameters, Kolmogorov-Smirnov statistics and its respectively 
p-values for the second data set. 



The good performance of our distribution is also confirmed by the Figure |2j Figures 
[3] and [4] show empirical versus fitted distribution functions of the models considered for 
the two data sets. Once again, we see that EPL distribution presents a good fit as the 
other distributions considered. 



7 Conclusion 

In this paper, a new lifetime distribution with decreasing failure rate, named exponential 
Poisson-Lindley (EPL) distribution, was introduced by compounding exponential and 





Figure 2: Plots of the fitted densities of the EPL distributions for the first (left) and 
second (right) data sets. 

Poisson-Lindley distributions. Several statistical properties of it were derived and dis- 
cussed, in particular: mean lifetime, moments, order statistics and Renyi entropy. More- 
over, maximum likelihood estimation and Fisher's information matrix were presented and 
discussed. Applications of the EPL distribution were carried out on two real data sets. 
Based on the results of these applications, it is shown that the EPL distribution provides 
a good competitor among the other well-known lifetime distributions. 
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Figure 3: Empirical versus fitted distribution functions for the first data set. 
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Figure 4: Empirical versus fitted distribution functions for the second data set. 



